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A number of deeply virtual exclusive experiments will allow us to access the Gener- 
alized Parton Distributions which are embedded in the complex amplitudes for such 
processes. The extraction from experiment is particularly challenging both because of 
the large number of kinematical variables and observables to be pinned down in each 
experimental analysis and because, at variance with inclusive experiments, the variables 
representing the quark momentum fraction appear integrated over in the physical am- 
plitudes and cannot be accessed directly. We present a strategy for the extraction from 
experiment that makes use of constraints from both elastic and inclusive scattering as 
well as information from lattice QCD results. 

1 Introduction 

Deeply virtual exclusive experiments such as ep — > e'^p' (Deeply Virtual Compton Scat- 
tering, DVCS) and ep e'Mp' (Deeply Virtual Meson Production, DVMP) can provide 
information on the partons' localization in space in addition to their longitudinal momentum 
fraction distribution [2 . The scattering amplitude for DVCS can be written within QCD 
factorization at leading order in l/Q as 

Q 

where is the four-momentum transfer squared in the hard collision, p{p') is the initial 
(final) proton momentum, t = is the four-momentum transfer squared between the initial 
and final protons, ^ = 2A+/(p-|-p') is the skewness variable, and 
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The GPD Fq{x,£^,t) (with Fq = Hq or Eq) is the soft part in the handbag diagram describ- 
ing this reaction, which is convoluted with the hard part, l/{x — £, + ie), and integrated 
over X, representing the partons longitudinal momentum fraction. Crossing symmetry is 
implemented by 

F,(±)(a;, t) = Fq{x, ^, t) T Fq{-x, t), (3) 

recalling that for PDFs, q{—x) = —<i{x) relates negative x to positive x antiquark probability. 

The goal of this contribution is to determine a physically motivated GPD parametriza- 
tion that satisfies known theoretical constraints from both elastic and inclusive scattering as 
well as information from lattice QCD results. The formulation that we present ultimately 
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aims at: "devising a form combining essential dynamical elements with a flexible model that 
allows for a fully quantitative analysis constrained by the data" . It is a step forward in the 
direction of current global analyses in related sectors such as the extraction of Transverse 
Momentum Distributions (TMDs) from semi- inclusive experimental data [3, 4 . Differently 
from TMDs however the extraction of GPDs is expected to be more challenging due to 
the large number of kinematical variables and observables to be pinned down in each ex- 
perimental analysis. In addition, as illustrated in Eq.([2|, at variance with inclusive and 
semi-inclusive experiments, the variables representing the quark momentum fraction appear 
integrated over in the physical amplitudes and cannot be accessed directly. Other features 
also appear that are unique to GPDs, and should be taken care of in a physically motivated 
parametrization, such as evolution (see e.g. discussion in [5j), the interplay between 
flavor separation and crossing symmetry [SJEJIT], the ability to constrain a given GPD from 
the limited number of available measurements [5] , and finally the applicability of dispersion 
relations [6l[7]. 

Here we discuss results obtained using the phenomenologically constrained parametriza- 
tion from Refs. [51 [TU] . and we present further developments of the model. The parametriza- 
tion is valid for unpolarized GPDs in the valence sector, at intermediate values of the skew- 
ness. It is aimed at having the same kinematical coverage of Jefferson Lab 6 GeV and 12 
GeV experimental programs. 

Our formulation is based on a quark-diquark picture improved by a Regge-type contri- 
bution (for more details see [9l [10] : 



ff(X,C,i) = G(X,C,t)i?(X,C,t), (4) 

where G and R are respectively the quark-diquark, and Regge-type contributions, and we 
dropped the subscript q for simplicity. Note that the parametrization is given in terms 
of the longitudinal variables X = jp^ , and C — /S.^ /p^ , which are in a one to one 
correspondence with the alternative set x,^ (see [llj for a review). The parametrization is 
constructed at an initial scale Q^, whose value is determined by our fit to the inclusive and 
elastic scattering data. H can subsequently be evolved perturbatively [T^ . 

Experimental constraints can be found by using the fact that H obeys the following 
relations [9j: 

H{X,Q,Q) = q{X) (5) 

1 

^^i/(X,C,i) = F,{t), (6) 

-i+C 

q{X) being the (valence) quark distribution, and Fi the Dirac form factor. Clearly, these 
constraints do not affect the C, dependence, however additional conditions can be found by 
using ab initio lattice QCD calculations of the higher moments in X ([TT and references 
therein). The n = 1,2,3 moments for both the isovector and isoscalr combinations: Hu-d = 
Hu — Hfi, and Hu+d = + are in fact available (moments for the GPD E are also 
available but we will not use them in this contribution) . Constraints from lattice QCD were 
used for the first time in Ref.[10 . 

By imposing the above constraints it was established in Ref.[n] that the quark-diquark 
picture cannot reproduce data/current parametrizations on q{X) at low X. One can in fact 
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expect that the quark-diquark or spectator models where the mass of the spectator is kept 
at a fixed value does not have the right physical input to reproduce the power dependence 
in X of the data. The term R{X, (, t) takes care of this problem effectively. This in turn 
allows for a quantitative description of both the low X behavior, and of the t dependence 
in the whole range of t. It is important to stress that a correct treatment of the low X 
behavior is crucial in determining the t dependence of GPDs, even if considering the valence 
region, because it dominates the form factor sum rule, Eq.(|6]). The latter cannot be satisfied 
without introducing a Regge type behavior at low X. This feature is missing from diquark 
models. A more thorough study is under way where Regge behavior emerges from a spectral 
distribution in the proton's soft debris invariant mass, Mx, valid at large values of the center 
of mass energy squared, s [T5]. 

For the present form of our parametrization we treated separately the so-called DGLAP 
{X > C) and ERBL {X < () regions. 

The functional forms for Eq.Q in the DGLAP region are given by 

G{xx,t) ^Y^J'^^^D{x,h^)D{xx,Ky ^ ^ 

R{XX,t) = ), (7b) 



where 



l-X ^ l-X 



M -Mi - {kj_ A , 

1-C l-X ^ V 1-C J l-X' 



In Eq.(7alA; and k' are the initial and final quark momenta respectively, D{X, \s.±) = k^~m? , 
D{XX,'k'±) = k'^ — m^, = k_L — (1 — X)/{1 — ()A, m being the struck quark mass. 
4>{k'^ , A) is the nucleon-quark-diquark vertex function, where the diquark is either a scalar or 



an axial- vector (this allows us to separate out the d and u quarks components). In Eq.(7bl 
tmin=-e/{l-OM\ 

Notice that the hybrid parametrization does not make use of a "profile function" for 
the parton distributions, but the forward limit, H{X, 0, 0) = q{X), is enforced non trivially. 
This affords us the flexibility that is necessary to model the behavior at ^, t ^ 0. The 
parameters in the DGLAP region were fitted separately: the set {a,A, Mx} was first used 
to determine the C, — Q,t — Q behavior, i.e. fitted to q{x) in the valence region; next, the 
set {p, 13} was used to fit to the nucleon form factor data. All parameters are listed in 
[S]. We reiterate that Q^, the initial scale, is also determined by the fit in our approach. 
More importantly, although all parameters are determined in our approach at ^ = 0, the C 
dependence in the DGLAP region is obtained through the kinematics expressed in Eqs.Q. 
The values of the parameters obtained for R do not correspond directly to similar ones from 
Regge theory because of the multiplicative form in Eq.Q that is such that e.g. the slope at 
low X receives contributions from both the R and G functions. 

To extend a parametrization in the ERBL region one needs C-dependent constraints given 
e.g. by the higher moments of GPDs that are calculable in ah initio calculations in lattice 
QCD. Our study is dedicated to defining and using these constraints. 
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In Figjl] we present one sample of our 
preliminary results from a study performed 
using the most recently available lattice 
evaluations from Ref . [14] . The figure shows 
n = 3 points for both the isovector, — 
—H^ , and isoscalar, H^ + H^, contributions 
plotted as a function of X at fixed kinemat- 
ics. The points where obtained by using 
a mathematical technique first used in |15j 
by which the function H(X, C, t) is recon- 
structed in average {H = H) over a given set 
of points, X in the interval X £ [0, C,], using 
a set of normalized polynomials, in this case 
the Bernstein polynomials, as weight func- 
tions around X {( and t are kept fixed) . The 
Bernstein moments, or the weighted aver- 
ages of H, are then obtained as linear com- 
binations of Ai{C,t), A2{C,t), As^^t) which 
are obtained from the lattice calculations in 
[14] . Notice that the latter are moments in 
the ERBL region, obtained by subtracting 
the DGLAP region from the total moments 
given in Ti]. Furthermore, the lattice calcu- 
lations on GPD moments need to be chirally 
extrapolated. In [10] we used a phenomeno- 
logical approach that can be applied to all 
n = 1,2,3 moments. The approach makes 
use of the fact that all moments seem to show a dipole behavior displaying increasing values 
of the dipole masses, A(") with increasing n. Chiral extrapolations for the dipole masses 
are then implemented. In [10] we checked that results were consistent with both data and 
alternative methods in the case n = 1. Here we were able to check that our new results are 
consistent with the newly available chiral extrapolations for n = 2 from Ref. [TB] . Our consis- 
tency check is important for testing the reliability of our method for n = 3 for which current 
calculations in chiral perturbation theory are impracticable, while well tested phenomeno- 
logical approaches like the ones proposed can be used. The errors on the x-axis represent 
the dispersion in our reconstruction technique while the errors on the y-axis are propagated 
from the lattice calculation errors. No other source of errors is taken into account. From 
Fig|l]one can see that the isovector contribution is consistent with previous determinations. 
One can notice a much improved error bar with respect to our previous analysis in [10 due 
to the improved set of lattice results. Furthermore, it is now possible to calculate also the 
isoscalar contribution. 

The dispersion along the the x-axis could improve if higher n moments could be calcu- 
lated using different lattice QCD techniques that avoid operator mixings and the associated 
renormalization issues [17 . Nevertheless, with better constraints in hand, such as the ones 
displayed in Figjl] showing both reduced errors on the isovector part, and for the first time 
the isoscalar contribution, one can now consider a fully fiedged parametrization in the ERBL 
region. 
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Figure 1: The isovector, Hu — Hd (top), and 
isoscalar, iJ„ -I- (bottom) contributions to 
the generalized parton distribution, H plot- 
ted as a function of X at fixed xbj — Ci t, 
and Q^, in the ERBL region (X < ( = 0.36) 
extracted using the lattice QCD calculations 
from Ref. [13]. The errors on the variable X 
represent the dispersion from the Bernstein 
moments method [TU], while the errors on H 
refiect the error on the lattice results. 
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In conclusion, we presented first results on a physically motivated parametrization that 
takes into account a number of constraints on the GPD H in the valence quarks sector. 
Similar results were also obtained for the GPD E. Parameters along with our statistical 
analysis are made available. Our work represents a step forward in the much needed direction 
of Global Parametrizations of both TMDs and GPDs. More comparisons with Jefi'erson Lab 
data on both DVCS and DVMP including the chiral-odd GPD sector are on their way. 
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